Curie temperature for an Ising model on Archimedean lattices 
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are 1.25, 1.40, 1.45, 2.15 and 2.80 [J /kg] for (3, 12^), (4, 6, 12), (4, 8^), (3, 4, 6, 4) and (3^, 6) lattices, 
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I. INTRODUCTION 

Beauty of the Ising model (IM) [J, Q manifests in its 
simplicity. The system considered is a network of N in- 
teracting spins Si = ±1 which energy is 

e^-jY,s,s,, (1) 

where J is an exchange integral. We assume homoge- 
neous short range spin interactions, i.e. the summation 
in Eq. |^ goes over pairs of nearest neighbours. 

Positive sign of J > gives ferromagnetic interaction 
among spins. The minimisation of energy for tem- 
perature T — Q produces spin dynamics which may be 
described by deterministic cellular automaton with rule 

5,(t + l) = sign|^J^^,(0j , (2) 

where t denotes discrete time and summation goes over 
nearest neighbours of i-th spin. 

For finite temperature T > the deterministic rule 
(0) is replaced by a probabilistic cellular automaton with 
spin update rule SAt) Si{t + 1) described by Glauber 
|3 or Metropolis 3 dynamics. Then the phase transi- 
tion may be observed: below critical temperature T < Tc 
spontaneous magnetisation m = J2iLi ^i/^ / is ob- 
served while m = for T > Tc- 

The IM was already investigated 

• in many ways, including Monte Carlo simulation 
, series expansion 0, Q , mean- field approach Q 
or partition function technique 0, 

• and for many systems, for example: antiferromag- 
nets flQ. 11. 12Lfrustrated (K), 12, 13], disordered 
[T^ or diluted jisll n etworks on complex ^j, or 
shuffled lattices |l7l|. etc. [l^- 



TABLE I: AL and associated critical temperatures Tc- 
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lattice 


Tc [J/ks] 


Ref. 


3 


(3,12^) 


1.25 






(4,6, 12) 


1.40 






(4,8^) 


1.45 






(6«) 


1.52 


[22] 


4 


(3,4,6,4) 


2.15 






(4*) 


2/arcsinh 1 ^ 2.27 


[20, 21J 




(3,6,3,6) 


2.27 


m 


5 


(3^6) 


2.80 






(3^42) 


2/ In 2 fa 2.89 


[25] 




(3^4,3,4) 


2.93 


[25] 


6 


(3«) 


3.64 


[22] 



In this paper the critical temperature Tc is estimated 
for five two-dimensional lattices, basing on {m{T)) de- 
pendence, where (• • • ) denotes the time average. The 
Archimedean lattices are vertex transitive graphs that 
can be embedded in a plane such that every face is a 
regular polygon. Kepler showed that there is exactly 
eleven such graphs [13 ■ The names of the lattices are 
given according to the sizes of faces incident to a given 
vertex. The face sizes are listed in order, starting with 
a face such that the list is the smallest possible in the 
lexicographical order. In this way, the square lattice 
gets the name (4, 4, 4, 4), abbreviated to (4^), honeycomb 
is called {6'^) and Kagome is (3,6,3,6). Some results 
concerning IM on AL were already presented in Refs. 
mmiillillllli^, however in a literature known to 
us the Curie temperatures of several AL are still missing. 

Critical properties of these lattices were investigated 
in terms of site percolation ^2^ in Ref. ■ Topologies 
of all eleven AL are given there as well. 



II. RESULTS OF THE SIMULATIONS 
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We evaluate the Curie temperature Tc basing on the 
termal dependence of magnetisation (m). Investigated 
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FIG. 1: Dependence of the average magnetisation (m) on 
temperature T expressed in [J/fcs] units for (S**, 6), (3, 4, 6, 4), 
(4,6,12), (4,8^) and (4*) AL. The simulations are carried 
out for TV 6 ■ 10* spins during Mtcr = 2-10^ [MCS]. The 
magnetisation (m) is averaged over the last 10^ [MCS]. 



systems contain about N k, Q ■ 10'* spins, which deco- 
rate nodes of (S'*, 6)^ (3, 4, 6, 4), (4, 6, 12), (4, 8^) AL. The 
Glauber dynamics 3] is appHed and the simulation takes 
iViter = 2-105 Monte Carlo steps (MCS). One MCS is 
completed when all N spins are investigated (spin-by- 
spin in a type- writer order), wheter they should flip or 
not. The time average is performed over the last 10'°^ 
MCS for an evaluation of (m) . The results are presented 
in Fig. ^ Temperature for which spontaneous magneti- 
sation (to) vanishes is accepted to be an estimation of 
Tc- These estimations are shown in Tab. ^ 



where k = T /Tc 



In contrast to Galam-Mauger [23| semi-exact formula 
for Tc dependence on system dimensionality d and lattice 
coordination number z, we show that critical tempera- 
ture for IM differ slightly for several AL (where d = 2) 
with the same values of z. Similarly to the percolation 
phenomena 26] , also for IM the dimensionalit y d and the 
coordination number z are not sufficient [30l l3ll Is^ for 
determining the critical point Tc- 
1 




0.9 0.95 

= T/Tc 

FIG. 2: Dependence of the average magnetisation (m) 
on normalized dimensionlees temperature T/Tc for (3^,6), 
(3, 4, 6, 4), (4, 6, 12), (4, 8^) and (4*) AL. The same data as in 
Fig. in 



III. CONCLUSIONS 



In this paper the Curie temperatures for IM on all AL 
are collected. Among them, Tc for (3,122), (4^6,12), 
(4,8^), (3,4,6,4) and (3'^,6) AL are evaluated for the 
first time with the Monte Carlo simulation. 

For all investigated AL, the shape of m{T/Tc) curve 
(see Fig. is roughly the same as for the square lattice. 
In the latter case, an analytical expression |2q is known 



1to(k)1 



'cosh'^(2/«;) 
sinh^(2/«;) 



[sinh2(2/K) - 1], 



Acknowledgments 



K.M. is grateful to Krzysztof Kulakowski for many 
valuable and fruitful discussions. Calculations were car- 
ried out in ACK-CYFRONET-AGH. The machine time 
on HP Integrity Superdome is financed by the Polish Min- 
istry of Science and Information Technology under Grant 
No. MNiI/HP.I_SD/AGH/047/2004. 



[1] W. Lenz, Z. Phys. 21 (1920) 613; E. Ising, Z. Phys. 31 
(1925) 253. 

[2] D.W. Heermann, Computer Simulation Methods in The- 
oretieal Physics, (Springer Verlag, Berlin 1990); D. Stauf- 
fer, F.W. Hehl, N. Ito, V. Winkelmann, J.G. Zabolitzky, 
Computer Simulation and Computer Algebra, (Springer 
Verlag, Heidelberg 1993). 

[3] R.J. Glauber, J. Math. Phys. 4, 294 (1963). 

[4] N. Metropolis, A.W. Resenbluth, M.N. Resenbluth, A.H. 



Teller, E. Teller, J. Chem. Phys. 21, 1087 (1953). 
[5] J. Garcia, J.A. Gonzalo, Physica A 326, 464 (2003); F. 
Merz, J.T. Chalker, Phys. Rev. B 65, 054425 (2002); N. 
Ito, K. Hukushima, K. Ogawa, Y. Ozeki, J. Phys. Soc. 
Japan 69, 1931 (2000); H.-O. Heuer, Annual Reviews of 
Computational Physics IV, ed. D. Stauffer, (World Scien- 
tific, Singapore 1996), p. 267; H. Rieger, Annual Reviews 
of Computational Physics II, ed. D. Stauffer, (World 
Scientific, Singapore 1995), p. 295; M. Acharyya, B.K. 



3 



Chakrabarti, Annual Reviews of Computational Physics 
I, ed. D. StaufFer, (World Scientific, Singapore 1995), p. 
107. 

[6] P. Butera, M. Comi, Phys. Rev. B 65, 144431 (2002); 
P. Butera, M. Comi, Phys. Rev. B 56, 8212 (1997); A.J. 
Guttmann, I.G. Enting, J. Phys. A 26, 806 (1993); J. 
Oilman, C.J. Hamer, W. Zheng, J. Phys. A 24, 2863 
(1991). 

[7] J. Adler, Annual Reviews of Computational Physics IV, 
ed. D. Stauffer, (World Scientific, Singapour 1996) p. 241. 

[8] F. Jasch, H. Kleinert, J. Math. Phys. 42 (2001) 52; J. 
Berges, N. Tetradis, C. Wetterich, Phys. Rev. Lett. 77 
(1996) 873; W.L. Bragg, E.J. Wilhams, Proc. Roy. Soc. 
A 145, 699 (1934). 

[9] A.K. Hartmann, cond-mat/0410583 K. Malarz, M.S. 
Magdon-Maksymowicz, A.Z. Maksymowicz, B. Kawecka- 
Magiera, K. Kulakowski, Int. J. Mod. Phys. C 14, 689 
(2003); A.F. Ferdinand, M.E. Fisher, Phys. Rev. 185, 
832 (1969). 

[10] M.J. Krawczyk, K. Malarz, B. Kawecka-Magiera, A.Z. 
Maksymowicz, K. Kulakowski, co nd-mat/0412545.. — to 
appear in Phys. Rev. B, 72 (2005). 

[11] J. Richter, J. Schulenburg, A. Honecker, D. Schmal- 
fufi, Phys. Rev. B 70, 174454 (2004); S.J. Penney, V.K. 
Cumyn, D.D. Belts, Physica A 330, 507 (2003); P. 
Tomczak, J. Schulenburg, J. Richter, A.R. Ferchmin, J. 
Phys.: Condens. Matter 13, 3851 (2001); B. Kawecka- 
Magiera, K. Kulakowski, A.Z. Maksymowicz, J. Magn. 
Magn. Mater. 221, 391 (2000); P. Tomczak, J. Richter, 
Phys. Rev. B 59, 107 (1999). 

[12] B. Tadic, K. Malarz, K. Kulakowski, Phys. Rev. Lett. 
94, 137204 (2005). 

[13] W. Lebrecht, E.E. Vogel, J.F. Valdes, J. Alloys Com- 
pounds 369, 66 (2004); 

[14] B. Kawecka-Magiera, A.Z. Maksymowicz, M. Kowal, K. 
Kulakowski, Int. J. Mod. Phys. C 9, 685 (1998); B. 
Kawecka-Magiera, A. Maksymowicz, K. Kulakowski, A. 
Lenda, Acta Phys. Pol. A 91, 297 (1997). 

[15] W. Selke, et al.. Annual Reviews of Computational 
Physics I, ed. D. Stauffer, (World Scientific, Singapour 
1995) p. 17; Z. Neda, J. Phys. I 4, 175 (1994). 

[16] A. Aleksiejuk, J.A. Holyst, D. Stauffer, Physica A 310, 
260 (2002); G. Bianconi, Physica A 303, 166 (2002); A. 
P§kalski, Phys. Rev. E 64, 057104 (2001). 

[17] K. Malarz, Int. J. Mod. Phys. C 14, 561 (2003); A.O. 



Sousa, K. Malarz, S. Galam, cond-mat/0503702" — to 
appear in Int. J. Mod. Phys. C, 16 (2005). 

[18] F. Mancini, New perspectives on the Ising model — to 
appear in Eur. Phys. J. B, (2005). 

[19] J. Kepler, Harmonices Mundi, (Lincii 1619). 

[20] H.A. Kramers, G.H. Wannier, Phys. Rev. 60, 252 (1941); 
H.A. Kramers, G.H. Wannier, Phys. Rev. 60, 263 (1941). 

[21] L. Onsager, Phys. Rev. 65, 117 (1944). 

[22] J.M. Dixon, J.A. Tuszyhski, E.J. Carpenter, Physica A 
349, 487 (2005). 

[23] W.T. Lu, F.Y. Wu, Phys. Rev. E 71, 046120 (2005); J. 
Zheng, G. Sun, Phys. Rev. B 71, 052408 (2005); R.J. 
Baxter, T.C. Choy, J. Phys. A 21, 2143 (1988); H. Gia- 
comini, J. Phys. A 21, L31 (1988). 

[24] K.Y. Lin, W.N. Huang, Aus. J. Phys. 38, 227 (1985); 
K.Y. Lin, J.M. Fang, Phys. Lett. A 109, 121 (1985); K.Y. 
Lin, T.L. Chen, Chin. J. Phys. 25, 178 (1987); K.Y. Lin, 
C.H. Kao, T.L. Chen, Phys. Lett. A 121, 443 (1987); 
K.Y. Lin, J.L. Chen, J. Phys. A 20 5695 (1987); K.Y. 
Lin, J. Stat. Phys. 49, 269 (1987); K.Y. Lin, S.C. Wang, 
Phys. Letters A 128, 143 (1988); K.Y. Lin, Phys. Letters 
A 128, 35 (1988); K.Y. Lin, Chin. J. 26, 150 (1988); K.Y. 
Lin, Int. J. Mod. Phys. B 3, 1237 (1989); K.Y. Lin, J. 
Phys. A 22, 3435 (1989). 

[25] C.J. Thompson, M.J. Wardrop, J. Phys. A 7, L65 (1974). 

[26] D. Stauffer, A. Aharony, Introduction to Percolation The- 
ory, (Taylor and Francis, London 1994); M. Sahimi, Ap- 
plications of Percolation Theory, (Taylor and Francis, 
London 1994). 

[27] P.N. Suding, R.M. Ziff, Phys. Rev. E 60, 275 (1999). 

[28] C.N. Yang, Phys. Rev, 85, 808 (1952). 

[29] S. Galam, A. Manger, Physica A 235, 573 (1997); S. 

Galam, J. Appl. Phys. 87, 7040 (2000). 
[30] S. Galam, A. Manger, Phys. Rev. E 53, 2177 (1996); 

S. Galam, A. Mauger, Phys. Rev. E 56, 322 (1997); S. 

Galam, A. Mauger, Eur. Phys. J. B 1, 255 (1998). 
[31] S. C. van der Marck, Phys. Rev. E 55, 1228 (1997); S. 

C. van der Marck, Phys. Rev. E 55, 1514 (1997); S. C. 

van der Marck, Phys. Rev. E 55, 3732(E) (1997); S. C. 

van der Marck, Int. J. Mod. Phys. C 9, 529 (1998); F. 

Babalievski, Phys. Rev. E 59, 1278 (1999). 
[32] S. Galam, K. Malarz, cond-mat/0408359 — to appear 

in Phys. Rev. E, 72 (2005); K. Malarz, S. Galam, Phys. 

Rev. E 71, 016125 (2005). 



